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Abstract 

We report exact results on the enumeration of close-packed dimers on a finite kagome 
lattice with general asymmetric dimer weights under periodic and cylindrical bound- 
ary conditions. For symmetric dimer weights, the resulting dimer generating func- 
tions reduce to very simple expressions, and we show how the simple expressions 
can be obtained from the consideration of a spin-variable mapping. 
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1 Introduction 



A central problem in lattice statistics is the enumeration of close-packed dimers 
on lattices and graphs. The origin of the problem has a long history dating back 
to a 1937 paper by Fowler and Rushbrooke [1] in an attempt of enumerating 
the absorption of diatomic molecules on a surface. A breakthrough in dimer 
statistics has been the exact solution of the generating function for a finite 
square lattice of size M x N, where M and N are arbitrary, obtained by 
Kasteleyn [2] and by Temperley and Fisher [3J in 1961. 

In view of the role of finite-size solutions in the conformal field theory discov- 
ered by Blote et al. |3] in 1986, it has been of increasing importance to consider 
solutions of lattice models for various finite two-dimensional lattices. Thus, the 
dimer solution has been extended to cylindrical [5J and nonorientable [6] lat- 
tices. However, these lattices are variants of the square lattice which may not 
necessarily exhibit special lattice-dependence features. 
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(a) (b) 

Fig. 1. (a) An M x N lattice of 6MN sites with asymmetric dimer weights, (b) A 
Kasteleyn edge orientation of the kagome lattice. A unit cell is the region bounded 
by broken lines. 

In a recent paper [7J we have reported enumeration results of close-packed 
dimers on an infinite kagome lattice with symmetric dimer weights (activi- 
ties) . The solution turned out to assume a very simple expression. In this pa- 
per we extend the solution to finite lattices with general asymmetric weights 
and under two different boundary conditions. We find the solutions given by 
entirely different expressions. For symmetric weights, however, the solutions 
again reduce to simple expressions. We show how the simple expressions can 
be deduced quite directly from the consideration of a spin-variable mapping. 



2 Finite lattices 



Consider a kagome lattice with asymmetric dimer weights x, y, z around up- 
pointing triangles and x', y\ z' around down-pointing triangles as shown in 
Fig. [m^a). We consider a lattice of M x unit cells having a total of 6MN 
lattice sites; the case of M = 1, = 3 is shown in Fig. [I](a). A unit cell of 
the lattice contains 6 sites numbered 1,2, ...,6 as indicated in Fig. [It^b). The 
Kasteleyn edge orientations adopted in [7] is also shown in Fig. [T](b). 

The method of Pfaffians [2] concerns with the evaluation of a 6MN x 6MN 
antisymmetric Kasteleyn matrix A written down according to edge weights 
and orientations (under specific boundary conditions) which can be read off 
from Figs[T](a) and[T]^b), and by adopting the prescription 

{+Wj^ orientation from i to ?' 
(1) 
—Wij orientation from j to i 

where Wij is the weight of edge ij. The dimer generating function is then 
given by the square root of the determinant of the matrix A. We consider the 
kagome lattice under two different boundary conditions. 
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2. 1 The periodic boundary condition 



First we consider the periodic (toroidal) boundary condition (PBC) for which 
the lattice is periodic in both the horizontal and vertical directions. Kasteleyn 
[2] has shown that under the PBC the dimer generating function ZpBc is a 
linear combination of four Pfaffians Pf|y4j|, z = 1, 2, 3, 4, 



^PBC 



- Pf |Ai| + Pfl^al + Pf IA3I + Pf|A4| 



Up to signs yet to be determined, the Pfaffians are the square root of the 
determinants specified by the Kasteleyn orientation of lattice edges with, or 
without, the reversal of arrows on edges connecting two opposite boundaries. 
A perusal of Fig. [T](b) and the use of the prescription ^ lead to the four 
QMN X QMN Kasteleyn matrices. 
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Here, the superscripts T denote transpose, ® is direct product, hi is the 
M X M identity matrix, and H]\f, T/v are the N x N matrices 



^ 1 ■ ■ ■ 0^ 
1 ■■■ 



H 



n 



■■■ 1 
-1 •■■ 



-N 



1 

1 



1 



^N 



1 

1 







(4) 



where Fn is to be used later in (ITSll . and 
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The determinant of a matrix is equal to the product of its eigenvalues. To 
determine eigenvalues of ^1,2,3,4, we first block-diagonalize the 4 matrices by 
appropriate Fourier transforms. Since and commute, they can be simul- 
taneously diagonahzed and replaced by respective eigenvalues e*^" and e~'^", 
where 6'„ = 27rn/N, n^O,...,N -1. 

Similarly, Tm and can be simultaneously diagonahzed and replaced by 
eigenvalues e*'^'" and e""^"", where 0^ = 27rm/M, m = 0, . . . , M — 1. 

Likewise, and Hj^ commute and they can be simultaneously diagonahzed 
and replaced by eigenvalues e*"" and e~*"", where a„ = (2n + l)7i/N. n = 
0, . . . ,N — 1; Hm and Hjj can be simultaneously diagonahzed and replaced 
by eigenvalues e*^"" and e"'^"", where = (2m + l)7r/M, m = 0, . . . , M - 1. 
Then we find 

det|Ai|= n n det \A{9n,(pm)\, det|A2|= [] 11 det|AK,0^)|, 

n=0 m=0 n=0 m=0 

N-lM-1 N-lM-1 

det|^|= n n det|A(e„,/3j|, det | A4I = [] 11 det |A(q;„, 

n=0 m=0 n=0 m=0 

(6) 

where the 6x6 matrix A{6, (p) is anti-hermitian and is given by 
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This yields 

detA{e, 

with 



= 2A + 2D cos((/)) + 2E cos(2^ + (p) + 4A2 sin^ ^ 

A = {xy'z + x'yz'Y + {x'y'z + xyz'Y 
D = —{xy'z — x'yz'Y 
E = {x'y'z — xyz'Y 
^2 = {x'yz — xy'z'Y. 
The desired generating function is now obtained by substituting either Pf|v4j 



/ det \Ai\ or — y det | Ai \ into ([2]) . In the present case the signs in front of the 
square roots can be determined by considering the case of M = iV = 1. By 
exphcit enumeration we have 



^PBC 



2{x + x'){yz' + y'z), 



M = N = 1. 



(9) 



It is readily verified that the expression is reproduced by ([2]) if all 4 terms 
in (E]) are positive. Thus, we are led to the final expression 



^PBC 



det I All 
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(10) 



where det lA,- 



1, 4, are given by 



For symmetric dimer weights x' = x,y' = y, z' = z, we have D = E = A2 = 0, 
det 1^1^2,3,4! = {4:xyz)^'^^ , and the simple result 



ZpBC = 2 ■ {4:xyz) 



MN 



x' = X, y' = y, z' = z. 



'Ill 



We shall see in Sec. [3] that this simple result can be understood and deduced 
directly using a spin- variable mapping. 
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In the case of an infinite lattice, ([2]) leads to the per-dimer free energy 
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+2E cos(e + 0) + 4A2 sin^ e 



(12) 



This free energy is independent of the boundary condition. The free energy 
( TT^ can also be deduced using the vertex- model approach introduced in 0, 
details of which are straightforward and will not be given. 

For symmetric dimer weights x' = x,y' = y, z' = z, (IT^ reduces further to 

f = ^ln{4xyz). (13) 

This result for an infinite lattice was first reported in |i8j with the full derivation 
given in [7j. The exact per-dimer entropy s = |ln2 obtained from (fT3l) at 
X = y = z = 1 has been cited earlier by Phares and Wunderlich from [9] and 
by Elser [10] from different considerations. 



2.2 The cylindrical boundary condition 



Consider next the cylindrical boundary condition (CBC) for which the lattice 
of M X unit cells is periodic in the horizontal direction. The Kasteleyn 
orientation is achieved by reversing the orientations of the 4M — 1 edges con- 
necting unit cells in the A^th column to those in the first column. This gives 
the dimer generating function as a single Pfaffian 

ZcBC = v/det|AcBc|, (14) 
where Aqbc is the 6MN x 6MN matrix 

^CBC =Oo,o ® hi (S) In + 0-1,0 ® ® Hn - «i^o ® hd ® Hn + '^0,1 ® Fm ® In 
- ao 1 ® Fm ®In + ® ^ Hn - 1 ® ® HJ^. 

(15) 

Here, matrices a are those in ([5]) and Fn has been given in (jl]). Again, the 
matrix Aqbc is block-diagonalized by replacing Hn and Hj^ by their respective 
eigenvalues. This leads to 

Af-l 

det lAcBcl = n det |SM(^n)|, e,, = (2n + l)7r/Ar (16) 

n=0 
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where Bm{0) is the 6M x 6M Kasteleyn matrix 



Bm{0) 



B®Im 

(be 





B_ B B+ 
B^ B 





















Fm + 5_ 

■ 

■ 

■ 








B B+ 
B_ B B^ 
B_ B 



(17) 



B = ao,o + e'^ai,o - e 
the 6x6 matrices 



1,0' B+ 



ao,i + e^^CLi,!, and B^ 



^0,1 



1,1 



are 



B 



y 




V 



z' 




-X 



ze 



xe 
-id 



-id 








-y 

x' + xe* 


-y 








-ze"- 

y 









-z' 




y X — xe 







-y' 

—x' + xe^^ 




B, 



^0 


ze-'^ 


-y\ 






































yO 





J 



:i8) 



( 0^ 




-ze^^ 

\ y oy 



The matrix Bm{0) is of a form of that occurring in the evaluation of an 
Ising partition function under the cyhndrical boundary condition [5|6] . and 
the determinant det \Bm(0)\ can be evaluated as follows: 



Let B^'^ denote the 6x6 matrix B with row i and column j removed. By 
Laplacian expansion and the use of a lemma established in [6] , the determinant 
of Bm can be expanded as 
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det I 5m I =det\B\- det \Bm-i \ + det ■ det \BI:j_^\ 

det \B^-^\ ■ det \B]^I^^\ + yze'^ det \B^'^\ ■ det \Blf_^\ 



+yze-'" det ■ det \B]:I_^\. 



(19) 



Similarly, the determinant of the matrix B]^} can be expanded as 



det l^l^^l = det ■ det \Bm~i\ + ^Met ■ det 

+1/2 det ■ det + yze^'^dei l^^'^'^'^l ■ det \B];}_^\ 

+yze-'^det ■ det \B]^}_iI (20) 

where B^'^'^'^ is the 6x6 matrix 5 with rows i and fc, and columns j and £ 
deleted. 



Write Bm = det |i?Af| and Cm = det|i?M|. Expansions (flQl) and fl20|) are 
recursion relations of B and C, 



— 0- Bm-1 + & Cm~i 
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(21) 



where 
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subject to the initial condition Bq = 1, Cq = 0. Explicitly using f[T8|) . fl22l 
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a =(a;2 + x'^){y''^z^ + y^z'^) + 4x2?/'22;'2 — 2{x'y' z + xyz'){xy'z + xyz) cos6 
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The recursion relation (|2T|) can be solved by introducing generating functions 

oo oo 

B{t) = J2 ^Mt^'^ C(t) = CMt"". (24) 

M=0 M=0 
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The recursion relation (12T1) gives 



B{t) = l + t[aB{t) + bC{t)] 

C{t)=t[cB{t) + dC{t)], (25) 
where we have made use of the initial condition Bo = 1, Co = 0. 

Solving ([25]) for B{t) and C(t), we obtain 

1 — dt 1 — dt 

" l-{a + d)t+{ad-bc)t^ " (l-A+t)(l-A_t)' ^'^^^ 

where 



A± =(a + d)/2 ± J(a - rf)V4 + be 



A + 2A2 sin^ e ± V (^ + 2A2 sin^ _ £)2 _ ^2 _ 2/}^ cos^ 0. 



(27) 



Here, A,D,i?,A2 have been given in ([8]). Partial fraction and expand the 
right-hand side of ( l26l) . and compare the resulting expansion with (Ell), one 
obtains 

\M+1 \M \M 

Finally, by combining (IT^ and (1161) . we obtain the desired generating function 



^CBC 



N-l 

n^A/(^n), en = i2n+l)n/N. 

n=0 



For symmetric weights x' = x,y' = y, z' = z we have 

a = 8x^?/^2;^(l — COS0), b = Sixy"^ z"^ sin 9 
c = Six^yzsmO, d = 8x'^y'^z'^{l + cos9) 

X+ = IGx'^y'^z^, A_ = 0. 

We obtain Bm{9) = (IGx^y^z'^y^ sm'^{9/2) and hence 

N-l 

ZcBC = {AxyzY"'' n = 2'-^ ■ (Axyz)''''. (29) 

n=0 

Again, we shall see in Sec. [3] that this simple expression can be understood 
using a spin-variable mapping. 

In the case of an infinite lattice and since A4. > A_, we have from (!28|) Bm{9) — > 
\¥ . This gives the per-dimer free energy 
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Fig. 2. (a) An A4 = 3 kagome lattice. Sublattice 1 sites are denoted by open circles, 
(b) The spin configuration deduced from a typical dimer configuration. The lattice 
is decomposed into 4 strips or loops depending on the boundary condition (see text). 
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where A+(^) is given by (p71) . It is readily verified that (!30|) is identical to the 
free energy (fT2l) after carrying out the integration over in (fT2l) . 



3 A spin-variable mapping 



The exact enumeration results ( fTTj) and ( !29l) for a lattice of M x unit cells 
with symmetric dimer weights are strikingly simple, suggesting the possibility 
of a simple derivation. Indeed, Zeng and Elser [TT] and Misguich et al. [12] 
have introduced a pseudo-spin consideration of enumerating quantum states 
which can be transcribed to the present x = y = z = 1 case flSj. However, the 
pseudo-spin consideration was presented in contexts of spin 1/2 antiferromag- 
nets and quantum dimer models, and the application to the classical dimer 
problem with general weights x, z is not immediately obvious. A simpler 
formulation is very much needed. 

We elucidate the matter by describing an alternate spin- variable mapping valid 
for general x, z. First we note that the numbers of x, ?/, and z dimers are 
always fixed for finite lattices. This is due to the fact that the three principal 
axes do not intersect at common points. The kagome lattice has three sublat- 
tices as numbered in Fig. [2|^a). Denote the number of sites on sublattice i by 
Aj, i = 1,2, 3, and the number of x dimers by A^., etc. Then as a consequence 
of the fact the principal axes do not intersect at common points, we have the 
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relations 

Ny + N, = N,, N, + N^ = N2, N^ + Ny = Ns. 

This gives rise to A^^; = (A'"2 + N-^ — Ni)/2, etc., which are fixed numbers. The 
dimer generating function is therefore a single monomial of the form 

so we need only to compute the constant fl. In the case of Ni = N2 = N^, = J\f 
we are considering, this leads to the expression 

Z = n{xyz)^/\ (31) 

where we have A/" = 2MN for both the PBC and CBC boundary conditions. 

We next map dimer configurations on the lattice to spin configurations on one 
sublattice, say, 1. Consider the example of the lattice shown in Fig. [2]^a). Let 
the lattice consist A4 rows of sublattice 1 sites and A4 + 1 rows of equal 2 and 
3 sites with open boundaries in the vertical direction. The boundary condition 
in the horizontal direction can be either open or periodic. Denote the number 
of sublattice 1 sites in the mth row by rim, which must satisfy the sum rule 

Til + n2 + h riM = Ni = even, (32) 

since the lattice must admit dimer coverings. 

Assign spin variable 

= ±1, m = 1,...,M, n = l,...,nm 

to sublattice 1 sites, where m is the row number counting beginning from 
the top. Adopt the convention that a = +1 (—1) if the dimer covering the 
site also covers a site above (below) the row. For +1 (—1) sites we remove 
the two edges below (above) the site as well as the horizontal edge directly 
above (below) it as shown in Fig. [2](b). This procedure decomposes the lattice 
into strips (loops) for open (periodic) boundary conditions in the horizontal 
direction. 

Now every strip (or loop) must have an even number of sites to accommo- 
date one (or 2) dimer covering(s). This condition imposes constraints on spin 
configurations that can be realized by this mapping. To best describe the 
constraints it is convenient to define a row variable 



llamn, m = l,2,...,M. (33) 



n=l 



It is then readily verified that we must have 
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ri = (-1) 



(-1) 



m = 



2,3,...,A^, 



Tm = 1. 



(34) 



The constraint tm = 1 is automatically satisfied due to fl32l) and the fact that 



TM = (n)(rir2) ■ ■ ■ (tm-iTm) = (-1)™^ + -+ = (-1)^^ = 1. (35) 



We remark that the row variable fl33|) can also be used to analyze the dimer 
models in higher dimensions considered in [14] . 

We can now compute the constant Q. Beginning with an overall spin state 
degeneracy 2^^ of sublattice 1 sites, each constraint in flMl) reduces the spin 
states by a factor of 2. Since there are Ai — 1 such constraints, we have 



Note that there is an extra factor 2-^~^^ for periodic boundary conditions in 
the horizontal direction since each loop has 2 dimer coverings. Expression (1361) 
is a very general result independent of specific values of rim- 

For a lattice of M x unit cells with toroidal boundary conditions PBC 
considered in Sec. EH we have Ni = N2 = N3 = J\f = 2MN, M = 2M. 
Hence (l36ll gives Vl = 2^^ ■ 2~(-^~^) ■ 2^^ where as in (!36|) the second factor 
is due to — 1 constraints with the A^th constraint automatically satisfied, 
and the third factor is due to the 2-fold dimer coverings of each loop. This 
leads to ^pbc = 2 • {ixyz)^^ in agreement with (fTTj) . 

For a lattice of M x unit cells with cylindrical boundary condition CBC 
considered in Sec. E^l we again have Ni = N2 = N3 = = 2MN, M = 2M. 
However, the A^th row of sublattice 1 spins must be all +1 (Cf. Fig. [2](b) 
with the bottom row of sites removed) reducing the counting by a factor of 
2-^ and the number of rows by 1. Hence n = 2-^ ■ 2^i-(^-i) ■ 2^ and 
ZcBc = 2^"^ ■ (Axyz)^^^ in agreement with (|29l). 

We remark that our spin-variable mapping is akin to one used recently by 
Dhar and Chandra [H]. However, the Dhar-Chandra approach focuses on an 
infinite lattice by ignoring what happens on the boundary. Here, we treat the 
boundary effect rigorously and apply the mapping to finite lattices. 




(M 1)^ open boundaries 
2Ni+2^ horizontal periodic boundary condition 



(36) 
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Abstract 

In this paper we continue our consideration of closed-packed dimers on the kagome 
lattice. Using the Pfaffian approach we evaluate the correlation between dimers on 
two lattice edges. It is found that the correlation is extremely short-ranged in the 
case of symmetric dimers weights. Explicit expressions for the nonvanishing correla- 
tions are obtained in the interior of a large lattice. We also describe a Grassmannian 
functional integral approach, and use it to evaluate the dimer generating function 
and correlation functions. 

Key words: kagome lattice, close-packed dimers, correlation function, 

Grassmannian approach 

PACS: 05.50.-Fq, 04.20.Jb, 02.10.Ox 



1 Introduction 

in tlie preceeding paper pQ, hereafter referred to as I, we presented exact re- 
sults on the generating function for closed-packed dimers on a finite kagome 
lattice with asymmetric dimer weights. To further illustrate the usefulness of 
the Pfaffian method used in I, in this paper we extend the consideration to 
dimer-dimer correlations. For symmetric dimer weights we find that the cor- 
relation is extremely short-ranged, a property unique to the kagome lattice 
and previously reported by us in [2]. Here we derive explicit expressions of 
nonvanishing correlation functions for a large lattice. We also describe the for- 
mulation of a Grassmannian function integral approach, and use it to evaluate 
the dimer generating funciton and correlation functions. 
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2 The dimer-dimer correlation function 



The dimer-dimer correlation function measures the correlation between dimers 
on two lattice edges. This correlation is best described by introducing an edge 
occupation number 



riij 



1 if edge ij is occupied by a dimer 
if edge ij is empty. 



Likewise, the edge vacancy number is fiij = 1 — riij. The correlation function 
between two dimers covering edge ij in unit cell at ri and edge M in unit cell 
at r2 is defined by 

C{ij, ri; k£, = {riij^n '^H.ra) - ('^ij.ri) (^^M.ra) 
= (nij,ri ?^H,r2) ~ {nij,ri){nu,r2) 

where (■) denotes the configuration average. 

The second line in ([2]) is useful in computing the correlation function in the 
Pfaffian approach [3], since using it we need only to keep track of the dimer 
generating function with specific edge(s) missing as dictated by {n) or (nn). 

Let A be the antisymmetric Kasteleyn matrix derived from a Kasteleyn orien- 
tation, and let A' denote the antisymmetric matrix derived from A with edge 
ij, say in computing (njj), missing. Write 

Z = PfA, Z' = Pi A' = Pi[A + A] (3) 

where A is the matrix with zero elements everywhere except the ij element is 
—Aij and the ji element is —Aji. Then 

(fiij) = Z'/Z = PiA'/PiA 

and 

, ,2 detA' det[y4(/ + G'A)] , 

where G = A^^ is the Green's function matrix and I the identity matrix. 

In (jlj) we need only to keep those row(s) and column(s) in A and A~^ where 
elements of A are nonzero. This effectively reduces dimensions of matrices G 
and A to at most 4 x 4 in the computation of ([2]). Explicitly, we have 
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('^ii,ri)^ = det 



I2 + 



G{ri,r2)ii G'(ri,r2) 



/ 







-A- 



\G{ri,r2)ji G{ri,r2)jj) \-Aj, 



where is the n x n identity matrix, G and A are the 4x4 matrices 



(5) 



G- 



A 



^G'(ri,ri),;,; ^(ri.ri)^^ G(ri,r2)ifc G{ri,r2)ii'^ 

G{ri,ri) ji G{ri,ri) jj G{ri,r2)jk G{ri,r2)ji 

G{r2,ri)ki G{r2,ri)kj G{r2,r2)kk G{r2,r2)ke 

\G{r2,ri)ii G{r2,ri)ij G{r2,r2)ik G{r2,r2)u J 

( -Aij ^ 

-Au 
-Afc 



-A,, 



The formulation so far is very general applicable to any finite lattice with 
symmetric or asymmetric dimer weights and ri and V2 arbitrary. We now 
specialize to a large lattice with symmetric dimer weights. 

In the interior of a large lattice, the correlation depends only on the difference 
r = ri — r2 = {tx-, r^y}, so elements of G are given by 



G{\rxx,riy\\ {r2x,r2yY)i 



Jo {2tiY ^ 



For symmetric weights, the 6x6 inverse matrix A ^( 
by using Equ. (7) in I, yielding 



(6) 

in ([6]) is computed 



A 



-1/ 



where 



3x3 



xye *^ 
-yxe^^ 



^ —zx zye 



-ie 



1 / ^3X3 + Q 



Axyz 



xz 

-yze^' 




3x3 



R. 



3x3 



-R 



3x3 



p. 



3x3 



^3x3 / 



Q 



3x3 



-ie 



X {e" — e"' 
yx 



(7) 



v 



zxe 



-ie 



—xy —xze 
-yz 
zy 



ie 
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Fig. 1. Labeling of the twelve edges belonging to one unit cell. 








xzie 



1) 
1) 



zy{l - e-*(^+'^)) 



V 



zxie 



and the superscript f denotes Hermitian conjugation. 

Substituting ([7j) into ([6]), we find G^ri^, rij^}; {r2x,r2y})ij = if |ria;-r2x| > 2 
or \riy — r2y\ > 2, or, explicitly, 



Equation ([8]) says that the correlation function vanishes identically if the dis- 
tance between the two lattice edges under consideration is larger than two 
lattice spacing. This is a consequence of the simple form of the free energy 
given by Equ. (13) in I. As pointed out in [2|, the absence of the dimer-dimer 
correlation beyond a certain distance, also found in the Sutherland-Rokhsar- 
Kivelson state of a quantum dimer model [US], is a property unique to the 
kagome lattice. 

We can now compute all nonvanishing correlation functions by substituting 
([6]) and ([7]) into ([2]), including correlations between two dimer edges within 
a unit cell and in two neighboring unit cells. There are 12 edges belonging 
to a unit cell as numbered in Fig. [TJ Details of the computations, which are 
straightforward, will not be given. We tabulate the results in Tables [T}131 



3 Grassmannian approach 

In a series of papers in 1980 [6j Samuel pointed out and explored the relation 
between the classical close-packed dimer problem and the free Grassmannian 
lattice field theory. He showed that the dimer generating function can be 



c(zj, ri; H,r2) = 0, |ri-r2|>2. 



8 
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represented as a fermionic path integral, and that dimer correlation functions 
are correlators occurring in a free fermion theory. He applied this formalism 
to dimers on the square lattice. Subsequently, Fendley et al. [7j extended 
the consideration to the triangular lattice using a slightly modified fermionic 
action; they also elucidated the field theory aspect of the formalism. Here we 
follow the formalism of [7] and apply it to the kagome lattice. 

We begin with the Kasteleyn orientation of the lattice shown in Fig. 1 in 
I. Associate real Grassmannian variables r/ to lattice sites and consider the 
functional integral 



exp 



Vt), (9) 



where Aij is the matrix element of the 6MN x 6MN Kasteleyn matrix A 
given by Equ. (1) in I and Vr] = Hi drji. 

It is well-known that the integral (Q is the Pfaffian of the matrix A, a fact 
which can be seen by expanding the exponential and using the result that the 
integral 



if ii ■ ■ ■ isMN is a permutation of 1, ■ ■ ■ , 6MN, and vanishes otherwise. Thus, 
the path integral (Q is identical to the dimer generating function. 

For finite lattices with toroidal boundary conditions PBC considered in Sec- 
tion 2.1 in I, the dimer generating function is again expressed as a linear 
combination of four integrals corresponding to periodic boundary conditions 
with, or without, the reversal of boundary edge orientations as discussed in 
Section 2.1 in I. . In each case the exponent in (Q can be block-diagonalize 
by Fourier transforms of r]i. 

For reverse edge orientations in both directions, for example, the Fourier trans- 
form is 



where (n, m) labels the position of unit cells, a = 1, . . . , 6 the sites within a 
unit cell, and = {2p + 1)it/N, p = 0, . . . , N - 1, Pg = {2q + 1)it/M, q = 



0,...,M- 1. 



After carrying out the Fourier transform, the functional integral becomes 

Z= fexpl^ V*ap,P„a M(^P,Pq) Vap,P„b]Vr], (10) 

where A[a,(3) is the 6x6 matrix given in Equ. (7) in paper I and Vr] = 
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Up,q,adriap,f3„a- Further using relations 



,a ~ 'la 



VaN-i-p,PM-i-q,ay 



and 



a,b 



ba 



— A{aN-l-p, Pm-i-q) 



ba 



we obtain from ffTOj) 
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where Vr]* = Up,g,adv, 



p,q,a "''lap,Pq,a 



Vr]* 



Y[p,q,a '^''7ajv_i_p,/3Af-l-g,a 



(12) 

Vt]. Here, we have 



used anticommutation relations of Grassmannian variables and renamed dummy 
variables p, q, a, b. 

Next we use a well-known formula of Gaussian integrals of Grassmannian 
variables (see, e.g., (2.8) in the first reference in [6]) to write the product of 
([TOD and ([H as 



p,g p,<i 

where Vrjp^gVr]*^ = I\adriap,Pq,adr]^^^p^^a- This is identically the expression 
det 1 given by Equ. (6) in I derived there using the method of Pfaffians. 

The Grassmannian formalism is most useful in evaluating dimer-dimer cor- 
relation functions. To evaluate the correlation function ([2]), for example, we 
note that the dimer generating function can be written as 



E n 



<ij> 



where riij is defined in ([T]) and the summation is over all dimer coverings of 
the lattice. Thus, we have 



nij) = WijZ' 



.1 dZ 



dwij 



-1 



d^z 



dwkfdw 



WijWkeZ / r]ir]jr]kr]eexp{A)Vr]. 



where A = (1/2) J^ijVi^ijVjy we have assumed that the Kasteleyn ori- 
entation is from i to j and from /c to £ to get the correct sign for the Grass- 
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mannian integrals. These are precisely expressions of correlators occurring in 
a Gaussian theory and can be immediately written down. This gives 

where A~^ is the inverse matrix of A with elements given by Equ. (1) in I. 
Here, we have used the Wick's theorem in expanding the four-point correlator. 

Finally, the dimer correlation function ([2]) is given simply by the expression 

c{tj; M) = w,,wu{-AllA-,l + A-.^A^]), (13) 

where A~^ = G for a large lattice has been given in ([H]). 

It can be verified that ffT^ gives rise to the same results as tabulated in 
Tables [T}131 Correlation functions of three or more dimers can be derived in a 
similar fashion. The advantage of using the Grassmannian method here is that 
it does not invoke products of large matrices needed in the Pfaffian approach. 
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Table 1 



Correlation function ([2]) between two dimers in the same unit cell, a = 1/16, b = 
-1/16, s = 3/16. 
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(r,,r,) = (0,1) 
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Table 2 



Correlation function ([2]) between two dimers in unit cells (0,0) and (0,1), a = 
1/16,6 = —1/16. Row indices 1,2, ...,12 label edges in unit cell (0,0); column 
indices label edges in unit cell (0, 1) (Cf. Fig. [T]). 
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(r,,r,) = (1,0) 
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Table 3 



Correlation function ([2]) between two dimers in unit cells (0,0) and (1,0), a = 
1/16,6 = —1/16. Row indices 1,2,..., 12) label edges in unit cell (0,0); column 
indices label edges in the unit cell (1,0) (Cf. Fig. [T]). 
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